A Critical "Dimension" in a Shell Model for Turbulence. 
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We investigate the GOY shell model within the scenario of a critical dimension in fully devel- 
oped turbulence. By changing the conserved quantities, one can continuously vary an "effective 
dimension" between d — 2 and d = 3. We identify a critical point between these two situations 
where the flux of energy changes sign and the helicity flux diverges. Close to the critical point the 
energy spectrum exhibits a turbulent scaling regime followed by a plateau of thermal equilibrium. 
We identify scaling laws and perform a rescaling argument to derive a relation between the critical 
exponents. We further discuss the distribution function of the energy flux. 



Many theoretical and experimental results for fully de- 
veloped turbulence have been offered over the last decade. 
A new approach has been presented by Yakhot ^ in 
which the method of generating functions by Polyakov 
is generalized to the Navier-Stokes equations. Apply- 
ing a renormalization group procedure results in an 
estimate of a critical dimension for turbulence, around 
dc ~ 2.5, thus following the foot steps of an original idea 
by Frisch and Fournier Q but correcting the actual value 
of the dimension. The physical idea behind the existence 
of a critical dimension is related to the well known fact 
that the energy cascade in three dimensional turbulence 
is "forward" (in fc-space) going from large to small scales 
whereas for two dimensional turbulence it is backward, 
from small to large scales. This leads to the identification 
of a critical dimension between two and three at which 
the flux of energy changes its sign, and the amplitude of 
the field turns into a peak where there is no flux neither 
forward nor backward, fn ref. |Q] the theory is expanded 
around this critical point in terms of a ratio between two 
time scales. However, it is not possible to investigate 
the physical behavior in a non-integer dimension directly, 
neither experimentally nor numerically. In this letter we 
therefore propose to study this type of criticality in a 
shell model for turbulence In particular we focus on 
the GOY model [^-^ which exhibits well known conser- 
vation laws: in the 3-d version energy and helicity are 
conserved; in the 2-d version energy and enstrophy are 
conserved. It is possible to continuously vary the effective 
dimension of the model by changing the second conserved 
quantity from a helicity to an enstrophy quantity. As the 
energy is always conserved, we can study the energy flux 
directly as a function of the variation in the second con- 
served quantity and we identify a critical point, where the 
flux changes sign. Indeed the second conserved quantity 
is non-physical at this point as expected. Nevertheless we 
are able numerically to extract a series of new properties 
of the spectrum and the PDF around this critical point. 
A similar observation of a change of sign in the energy 
flux as a function of a parameter was already made in 
a different shell model by Bell and Nelkin In their 
model the dynamics is not intermittent and the proper- 
ties of the model are thus quite different from the GOY 
model. 



Our starting point is the approximative approach to 
turbulence made by discretizing the wave number space 
by exponentially separating "shells" In this respect, 
we apply the GOY model which has been successful 
in giving results for intermittency corrections in agree- 
ment with experiments H] (for other results on the GOY 
model, see The starting point is a set of wave 

numbers fc„ = fco2" and an associated complex ampli- 
tude Un of the velocity field. Each amplitude interacts 
with nearest and next-nearest neighboring shells and the 
corresponding set of coupled ODE's takes the form: 
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with n = 1,---N, kn = r" kg {r = 2), and boundary 
conditions bi — bpj — ci — C2 — apj^i = gn = 0. 
The values of the coupling constants are fixed by impos- 
ing conserved quantities. By conserving the total energy 
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when / 



— 0, we obtain the constraints 
On + bn+i + c„4.2 = 0. The time scale is fixed by the con- 
dition a„ ~ 1 leaving free the parameter S by defining 
the coupling constants as 
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The model also possesses a second conserved quantity of 
the form 



(3) 



which leads to a relation between a and 6: 2" — l/{6—l). 
For 6 < 1 this relation requires complex values of a, 
with Q{a) — n/ln 2. In 3d turbulence helicity H = 
/(V X u(x))-u(x)dx is conserved, which in terms of shell 
variables takes the form 
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when the values of parameters are S — ^ and K(q;) = 1 . 
In 2d turbulence on the other hand enstrophy — 
/ |V X u(x)pc?x is conserved and this corresponds to the 
parameters 6 = |, a = 2 which on the shells takes the 
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form ^ = J2 Kl'^nl'^ ■ Note that energy is conserved for 
any value of 6. This gives us the possibiUty to continu- 
ously vary the effective dimension (and thus the general- 
ized helicity/enstrophy (^) by varying the parameter S 
between 6=^ (3d) and (5 = | (2d). 



/). According to Eq. (^), the generalized helicity Q di- 
verges at 5 = 1 (i.e. a oo) and this could be a reason 
for the inhibition of the energy transfer [Q. Further- 
more, the rate of injected generalized helicity, given by 
E„(-l)"^«3^(/nO, diverges at this point 
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FIG. 1. a): Average energy flux versus the wavenum- 
ber k with N ^ 25,u = 10''^°, Uf = 2 for 5 = 0.5(x), and 
e = 2 • 10~'^{+), 2 • 10"^(*), . . ., 2 • 10"*(D). Note that as e 
decreases the inertial range shrinks, b); Inverse energy flux 
with A = 33, = lO"^*^ for e = -0.125. The forcing term is 
on shell 15 and of the form /15 = (1 -I- i) * 10~*/ui5. A large 
scale viscosity is now applied, see the text. 



The critical point is identified by looking at the energy 
flux through each shell which is given by [p3 



dt 

-/to ( fc„u„u„+i(w„+2 + ^—^r-^Un-i)) ) ) , (5) 



where only the contributions of the nonlinear terms are 
considered in the time-rate-of-change of the cumulative 
energy. From Eq. ^ we see that the last term vanishes 
as (5 ^ 1 causing a depletion in the energy transfer . 
This is observed in the numerical simulations shown in 
Fig. 0a, where the inertial range of the flux shrinks as 
e = 1 — S Q. Note, that we in Fig. la apply a forcing 
on the form /„^. = (1 -|- i) * 10~^/m*^ in order to ensure a 
constant input of the energy (and thus a constant flux). 
We reach similar conclusions when a constant determin- 
istic forcing is applied. Moving above 6=1, the energy 
flux reverses going instead from small to large scales, see 
Fig. 0b. Therefore the point S = Sc = 1 defines a criti- 
cal point where the energy flux for finite value of energy 
input / discontinuously jumps from positive to negative 
values (the jump diminishes with the forcing amplitude 
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FIG. 2. a): The spectrum < [m| > versus k for the 3d 
case S = 0.5(0) and for e = 2 • 10"''(+), . . ., e = 2 • 10~"'(*). 
Note the flat part of the spectrum developing as e — > 0. b) : A 
rescaling plot of the spectra in Fig. (2) by < \ u\ >*(^)'''^* 
versus fc*(^)"°-^ (e = 10"''(o) and e = 10"'^ (A)), where 
e,. = 2 • 10"^(*). 

Let us now turn to the spectra for 6 < 1. As (5 is in- 
creased from (5 = 0.5 and e ^ 0, one observes a build up 
of a shoulder leading to a plateau in the spectrum at large 
values of fc„. This is shown in Fig.^j, where the parameter 
e is varied one decade for each spectrum (note that in our 
spectra we plot < |m„| > vs. fc„ which provides similar 
information as plotting < |m„P >). We identified various 
scaling laws associated with the spectra of Fig.||. First 
of all the dissipative cut-off, fc^, moves in a systematic 
fashion as a function of e. We find the following scal- 
ing law: kd ~ e""* , with ad — 0.3. The plateau in the 
spectra (equipartition of energy among the shells) may 
be interpreted as a thermal equilibrium which overcomes 
the turbulent regime when the forward transfer of energy 
is reduced. We found that the level of the plateau scales 
with £ as (I Un \)pi ~ e~"p' , where api ~ 0.28, thus finally 
turning into a diverging amplitude around the forcing 
scale. The turbulent, cascading, part of the spectrum 
varies like (|u„|) k""" , as — 0.4, taking into account 
corrections due to intermittency Finally, the critical 
wavenumber kc, at which the spectrum crosses over from 
turbulent behavior to thermal equilibrium, also moves 
with e, possibly like kc ^ e""- To determine kc, we bal- 
ance the contributions from the two regimes 



{\Un\) {\Un Dpi 



(6) 



2 



and obtain the following scaling law for fee 



^% Qic = api/as ~ 0.7 



(7) 



showing that the scaling exponents are not all indepen- 
dent pq ] . This result can be verified by a simple rescaling 
of data. Let us assume that {\un\) / {\un\) pi is a function 
of k/kc alone, i.e. 



{\Un\)pl kc 
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where f{x) is such that f{x) ^ a;~"=, x << 1 and 
f{x) ~ const, X » 1. Then a data collapse is obtained 
by plotting {\un\) /e^°'p' versus k/e"". A good rescaling 
plot is obtained, see Fig||b, when the estimated value 
ac = 0.7 is used. Notice that the collapse of data does 
not apply to the dissipative range, since kd and kc scale 
differently with e. 

Since there is no transfer of energy at the critical point, 
the non-linear terms will not play any role and the equa- 
tions will only include the dissipation and forcing terms. 
This can be made quantitative by the fact that the state 
with a peak at the forcing scale 
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is a fixed point of the equations (0) which at e = is 
marginally stable p^ . Indeed we find numerically at 
e = that by starting with the fixed point (||) the peak 
stays at the forcing scale and the amplitudes remain zero 
above but become non-zero, although small, below the 
forcing scale. On the contrary, for e — > 0^, the peak 
is unstable and the energy is soon redistributed to the 
neighbouring shells. The existence of a sharp transition 
into the critical point was already indicated by a calcu- 
lation of the maximal Lyapunov exponent which drops 
sharply as e ^ In order to study the behavior of 
the spectrum around the critical point in details we have 
two parameters to vary, the "dimension" parameter, (5, 
and the viscosity, v. Fig. ^ shows a series of spectra for 
e — 0.002 varying v. Again one observes the shoulder 
at large k. As expected, the shoulder moves to higher 
k when the viscosity decreases. The wavevector for the 
dissipative cut off A:_d moves in the Kolmogorov fashion 
kd v~^/^. This leads us again to perform a standard fi- 
nite size rescaling plot, rescaling the k axis by {v/viof'^'^ 
and the velocity axis by (z^/z^io)~"'^*, see Fig. ^j. 
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FIG. 3. 
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a): The spectrum < | u | > versus k for the stan- 



dard 3d case - 5 — 0.5, v 
&nAu = 10""'(A), 
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\a) - and for e = 2 • 10" 



,10 ^(+). b): A rescaling of the curves 



in a) by < I u I > *{u/vio) 
10-"'. 



k * {{u /vxof'^^ where 



Now consider the other side of the critical point, 
namely on the "two-dimensional" side for 5 > 1. The 
behavior of 2d shell models has been previously inves- 
tigated in several papers [p"8|-p0[. A kind of "coupled 
GOY model" [|l^ gives an inverse fiux of energy which 
is explained in terms of a mean diffusive drift in a sys- 
tem close to statistical equilibrium, and shell models 
for 2d do not seem to give an inverse energy cascade 
with the usual "5/3" spectrum. In order to extract the 
energy of the inverse cascade we need to add a large 
scale viscosity to Eqs. (|l|) of the type —ly'k^'^Un- In 
Fig. (^) we show the energy spectrum for the case 
S = 1.125, 10-^^, = 10"^ The two branches of 
statistical equilibrium, respectively energy and general- 
ized enstrophy (^) equilibrium are clearly visible (see pC| ] 
for details). 

Let us turn our attention to the probability density 
functions (PDF). It is well known that, in fully developed 
turbulence, the PDF at the largest scales (small k) typi- 
cally behaves like a Gaussian, slowly changing its form as 
one moves towards the small scales (large fc), turning into 
a shape where large events play an important role giving 
a kind of stretched exponential PDF. Such a non trivial 
behavior is related to the property of multiscaling of the 
structure functions. When a turbulent scaling regime is 
detectable as in Fig. ^, we observe that the intermit- 
tency corrections appear to persist, even close to e = 0. 
However, in the flat part of the spectrum of Fig. ||, the 
probability distribution becomes wider at larger k but 
such that the shape is simply rescaled onto a universal 
and almost symmetric curve (rescaling by the standard 
deviation), see Fig. ||. 
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FIG. 4. <\u\> versus k ior 5 = 1.125 (a = 3). The two 
branches of statistical equihbrium < | u | > ~ const (energy 
equihbrium) and < | u j > ~ k~^^^ (generalized enstrophy 
equilibrium) are clearly visible. 
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FIG. 5. Rescaling of the PDF obtained for shells 
in the "flat" part of the spectrum in Fig. 2 with 
f = 10"^°, e = 3 ■ 10"^ for the shells n = 4, 6, 8, 10, 12. The 
pdf's are in each case rescaled by the standard deviation. 



In this Letter we have presented results for the exis- 
tence of a critical point in the "GOY" shell model. Our 
main results can be summarized as follows. At this criti- 
cal point, which lies between three- and two-dimensional 
behavior, the energy flux changes its sign, going from a 
forward to a backward transfer. Approaching this point 
from the "three-dimensional" side, part of the spectrum 
becomes flat as an indication of a thermal equilibrium. 
The cross-over to the flat part is determined by balanc- 
ing the turbulent energy spectrum with a spectrum in 
thermal equilibrium. We identify the scaling behavior of 
the cross-over point and rescale the spectra accordingly. 
The PDF of the thermal equilibrium shows simple scaling 
invariance although the statistics is not Gaussian. For 
analytical understanding of these results, one rewrites 
the "GOY" equations in terms of a generating function 
technique thus obtaining a set of coupled ordinary 
differential equations [|l|. One can further map these 
equations onto a Fokker-Planck equation for the distri- 
bution of the exponentiated quantities. We will discuss 
this in a forthcoming publication. 
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3 +C{5). Since 



the turbulent spectrum is infared divergent the main con- 
tribution to the energy will be associated the forcing shell 
Uf. The equilibrium spectrum is on the other hand dom- 
inated by the dissipation scale kd so the total energy 

2 

is _E = — ^kn, + C{5)kd- As the dissipation scale 
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wiU vary like kd ~ {&c — 5)* then C{S) ^ {5c — 
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£(fc„) = — ^fc-i +0((<5c-<5)-A, fee ^ ((5c -<5)*. 

(Sc~S) 3 
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